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Abstract 
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is a sum of 2- and 3-body interactions with no undetermined parameters. 
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pion photoproduction on the deuteron is in agreement with experiment. 
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I. INTRODUCTION 



In recent times, there has been increasing interest in applying the method of chiral 
lagrangians, or chiral perturbation theory (xPT), to processes involving more than a 
single nucleon [1-8]. This interest is motivated by the desire to determine what aspects of 
nuclear physics can be understood on the basis of the chiral symmetry of QCD. In the 
limit of vanishing u and d quark masses the QCD lagrangian admits a global SU(2)l X 
SU(2)r chiral symmetry. The absence of parity doubling in the hadronic spectrum implies 
that this chiral symmetry is either anomalous, or spontaneously broken by the QCD 
vacuum. The existence of suitable Goldstone boson candidates, the pions, bears out the 
latter conjecture. A priori, the space of chiral symmetric operators involving the relevant 
degrees of freedom — nucleons and pions — is infinite. Moreover, since the interactions 
are strong, nothing is really learned insofar as QCD is concerned by restricting oneself 
to any arbitrary finite subset of chiral symmetric operators; e.g., the simplest operators 
involving the fewest number of fields, or renormalizable operators. Fortunately, as a 
consequence of non-linearly realized chiral symmetry, low-energy hadronic matrix elements 
involving pions are analytic in momenta. The parameters that appear at leading order in 
small momenta are well known, and so chiral symmetry leads to non-trivial predictions 
at energies near the threshold of a physical process — so-called low-energy theorems 
(LETs). The crucial fact that the u and d quark masses are not identically zero can 
be straightforwardly accounted for in perturbation theory. As with any approximation 
scheme, the fundamental importance of xPT is its ability to handle corrections to the 
leading order in a systematic way ||. This method has been applied with great success 



to the interactions of pions with a single nucleon [10] [11|. One might then wonder what 



generic features of nuclear physics can be deduced from this chiral symmetry of QCD. 

The main technical difficulty that arises when considering more than a single nucleon 
is that xPT necessarily breaks down, as is made clear by the appearance of shallow 
nuclear bound states [[J. This breakdown manifests itself via infrared singularities in 
Feynman diagrams evaluated in the static approximation. The problem is clear in the 
language of time-ordered perturbation theory. Evidently there are two types of energy 
denominators that can appear in a typical time-ordered graph. The first type arises from 
intermediate states which differ in energy from the initial and final states only by the 
emission or absorption of soft pions. These energy denominators are of the order of a 
small momentum or pion mass and are therefore consistent with the usual chiral power 
counting scheme. On the other hand, the second type of intermediate states differ in 
energy from the initial and final states by a nucleon 3-momentum, and therefore blow up 
in the static limit. Graphs of the first type are called irreducible. Following the tenets 
of scattering theory, one can modify the rules and use xPT to calculate an effective 
potential, which consists of the sum of all iV n -nucleon irreducible graphs [[I] 0. The 
S-matrix, which of course includes all reducible contributions, is then obtained through 
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iteration by solving a Lippmann-Schwinger equation. Several generic features of nuclear 
physics, including two- and few-body forces M M |J and isospin violation II, have been 
shown to arise naturally in this approach. 

Here we will apply chiral perturbation theory to a scattering process involving light 
nuclei. This program was initiated several years ago by Weinberg |3|], in a study of the 
pion-nucleus scattering lengths. The strategy is best described graphically. In figure 1 
we display the anatomy of a scattering matrix element. xPT generates the irreducible 
kernel, which is then sewn to the external nuclear wave functions (with N n = A). Nuclear 
wave functions, too, are calculable in x^T, and in fact exist for the deuteron Q. A 
completely consistent calculation can be carried out, and results are reported here. How- 
ever, this approach is in a sense counterproductive, since it is then unclear what aspect 
of chiral symmetry is being tested. In the spirit of X-PT, where experimental input is 
always welcome, one should make use of the most successful phenomenological potential. 
This allows one to test the relevance of chiral symmetry in determining the irreducible 
scattering kernel. In the case of the deuteron, we also use the well known Bonn potential 
wave function ||12|| , which although not respectful of chiral symmetry, is at least spiritually 
linked to QCD. 

In this paper, we consider pion photoproduction on light nuclei at threshold. We carry 
out a general power counting analysis of photoproduction, but specialize to neutral pion 
photoproduction. We calculate a general formula for the invariant threshold amplitude 
to third order in small momenta. This process is of course interesting in its own right; 
there has been no systematic calculation based on chiral symmetry. Furthermore, as we 
will see, this process is intimately related to threshold neutral pion photoproduction on 
a single nucleon, a process which has caused a great deal of theoretical and experimen- 
tal confusion. In contrast to pion-nucleus scattering at threshold the amplitude for 
neutral pion photoproduction on nuclei, to third order in small momenta, has no undeter- 
mined parameters. Unfortunately, although the same is true in the single nucleon sector, 
evidently the amplitude there converges slowly at best. Hence the single-scattering con- 
tribution must be treated as phenomenological input. With a reasonable choice for this 
contribution, our result for the deuteron electric dipole amplitude is in agreement with 
experiment. 

This paper is organized as follows. In section 2 we review the standard power counting 
formulas. In section 3 we present the heavy fermion effective lagrangian to the order 
relevant to our calculation. Section 4 consists of a power counting analysis of pion photo- 
production on a nucleus composed of A nucleons. Here we give a general formula for the 
neutral pion photoproduction amplitude . We specialize to neutral pion photoproduction 
on the deuteron in section 5. Finally, we summarize and conclude. Several appendices 
are included for pedagogical purpose. 
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II. POWER COUNTING 



In the absence of exact solutions to quantum mechanical equations of motion, as in 
QCD, systematic statements are possible only when a small dimensionless expansion 
parameter is identified. One class of dimensionless parameters consists of pure numbers, 
such as a coupling constant associated with a renormalizable interaction, or the inverse 
of the dimensionality of a group, as in the large- M c limit. A second class of dimensionless 
parameters consist of ratios of dimensional quantities. This class arises naturally via 
broken symmetries, and is the case of interest in this paper. 

Spontaneously broken continuous symmetries give rise to massless Goldstone modes. 
These modes do not propagate in the vacuum, and therefore only couple derivatively. 
Hence at energies small relative to the characteristic symmetry breaking scale, the in- 
teractions of Goldstone bosons admit a power series in momenta. In QCD the small 
parameter is Q/A x , where Q is a characteristic momentum, and A x is the scale of chi- 
ral symmetry breaking — of order the masses of the lowest lying resonances. The effects 
of non-zero quark masses give the pion a mass. Since M n /A x is small, we treat Q as 
representing a small momentum or a pion mass. A generic matrix element involving the 
interaction of any number of pions and nucleons can then be written in the form 

M = Q v HQ/^)i (1) 

where \i is a renormalization scale, and v is a counting index. It is straightforward to 
arrive at a general formula for v by considering the momentum space structure of generic 
Feynman rules. In this way one finds 

v = 4L-I n -2I n + + 4 - (2) 

where L is the number of loops, 7 n>7r is the number of internal fermion, pion lines, V{ is 
the number of vertices of type i, di is the number of derivatives or powers of M n which 
contribute to an interaction of type i, and C is the number of separately connected pieces. 
One can also make use of the graph theoretic identities: 

L= I n + 4 - YXi + 1 (3) 

i 

2/ n + E n = Y^ViTH, (4) 

i 

where E n is the number of external nucleon lines, and is the number of fermion fields 
involved in an interaction of type i. Here we are interested in processes with the same 
number of nucleon lines in the initial and final state, and so defining N n = E n /2, we 
obtain the master formula 
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v = 4 - N n - 2C + 2L + YfyAi 

i 

Aj = di + Ui/2 - 2. 



(5) 



This formula is important because chiral symmetry places a lower bound: Aj > 0. Hence 
the leading irreducible graphs are tree graphs (L = 0) with the maximum number C of 
seperately connected pieces, constructed from vertices with A, = 0. 

How is this analysis altered in the presence of an electromagnetic field? Photons couple 
via the electromagnetic field strength tensor and by minimal substitution. This has the 
simple effect of modifying the lower bound on Aj to Aj > —1. (And, of course, of 
introducing an expansion in the electromagnetic coupling.) 



III. BARYON X PT 

With the power counting scheme established, the next step is to construct the various 
interactions which contribute to matrix elements for a given value of v. The technology 



that goes into building an effective lagrangian is standard by now |I(J . Here we establish 
our conventions. The pion triplet is contained in a matrix field 

£ = exp(^— -), (6) 

which transforms under SU(2)l x SU(2)r as S —>■ LUR) . It is convenient to introduce 
the field £ = vE, with transformation property £ — > L^W = U^Ftf. This transformation 
property implicitly defines U. Out of £ one can construct 

V, = \[^{d, - ieA„Q)i + C{d, - ieA»Q)g] (7) 
^ = - ieA„Q)Z - - ieA^Q)?], (8) 

which transform as V„ -> UV^W + Udjfl and A M -> UA„W under SU(2) L x SU(2) R . 
It is convenient to assign the nucleon doublet A the transformation property A — > C/A. 
With these ingredients, one can construct the leading order effective lagrangians, 

4 2 2 = \f>{D^D^) + ~^ 2 M>(£ + St) (9) 
4v = iN(P~ m)N + g A N4 l5 N (10) 
A = ^(AIW) 2 , (11) 

where = <9 M + V^, T a is an arbitrary Hermitian operator, and the C a are undetermined 
coefficients. The pion covariant derivative is D^L = — ieA^[Q, £], where Q = (1 + 
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Tz)/2. The appearance of the nucleon mass m, both explicitly and through the time 
derivative acting on the nucleon field, implies the existence of a dimensionless quantity 
that is not small: m/A x ~ 1. This destroys the power counting. Fortunately, since 
the nucleon carries a quantum number — Baryon number — which is conserved by the 
strong interactions, one can maintain a consistent power counting framework by choosing 



a heavy-fermion basis in which the nucleon mass does not appear at leading order [13 



[TJj]. The nucleon momentum can be written as p M = mt> M + where is the nucleon 



four- velocity (v 2 = 1), and k 11 is the amount by which the nucleon momentum is off shell. 
We can then define a velocity dependent basis 

B v {x) = e imv - x N(x). (12) 

In the rest frame of the nucleon, w M = (1, 0, 0, 0). Although choosing a particular velocity 
breaks covariance, integrating over all velocity-dependent lagrangians restores covariance 
|13| . The velocity-dependent lagrangian in the new basis is 

£« = iB v {v ■ D)B V + 2g A B v {A ■ S V )B V , (13) 

where S% is the spin operator. Note that the nucleon mass term is no longer present in 
the lagrangian. The new effective expansion parameter is \k\/A x . The 1/m corrections 
enter via higher dimensional operators which affect the higher orders in xPT . At next to 
leading order one finds, 

= ^-Bv{~D 2 + (v ■ Df + 2ig A {v ■ A, S ■ D} 
2m 

- l -[S», S:][(l + K v )f+ + 1 -{k s - K v )trf+] + ...)B V , (14) 

where k v = n p — K n , k s = K p +K n , and /+ = e^Q^+^Q^^F^. F^ u is the electromagnetic 
field strength tensor. The dots signify that we have included only the interactions at this 
order which are relevant to our calculation. Note, finally, that the A isobar can be 
introduced in an analogous way. 



IV. PION PHOTOPRODUCTION ON LIGHT NUCLEI 

Now we can put our technology to use. Consider a process with N n = A nucleons in 
both the initial and the final state, and a single photon and a single pion in the initial 
and final state, respectively. Limiting ourselves to lowest order in the electromagnetic 
coupling, we can order the chiral expansion of the irreducible diagrams by way of the 
counting index v: 

v = 3 - 3A 
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At leading order in small momenta, the matrix element is given by tree graphs with the 
maximum number of separately connected pieces (L — 0, C — A) constructed out of one 
interaction with A, = — 1, and interactions with Aj = 0. For example, at this order one 
has the the Kroll-Ruderman term (figure 2a). 

v = 4 - 3A 



The first corrections to the leading terms are still tree graphs with the maximum number 
of separately connected pieces (L = 0, C = A), but have one vertex of higher index: i) 
either one vertex with Aj = — 1 and one with Aj = 1; ii) or vertices with Aj = 0, the one 
involving the photon field being a 1/m correction (figure 2b). 

v = 5 - 3A 



There are four classes of corrections at this order: 

(i) One loop graphs (L — 1) with interactions with Aj = —1,0, and C = A (figure 2c). 

(ii) Counterterm graphs with L = 0, Aj = 1 and C = A (figure 2d). In the case of 
neutral pion photoproduction, these graphs vanish and only finite loop graphs remain. 

(iii) C = A tree (L = 0) graphs with i) Aj = —1,2 interactions; or ii) Aj = 0,1 
interactions (figure 2e) f\ (Some of these are 1/m 2 corrections and are proportional to the 
nucleon magnetic moments.) 

(iv) Finally, there are tree graphs (L = 0) with one less than the maximum number of 
separately connected pieces (C = A — 1), and interactions with Aj = —1,0. These graphs 
fall into two separate classes. There are the 3-body graphs like the Feynman graph of 
figures 2f and the time-ordered graph of 2g. For A > 2, there are also disconnected 2-body 
interactions as in figure 2h. 

In principle, all of these graphs contribute to pion photoproduction on nuclei. However, 
some generic simplifications arise when one sews in the nuclear wavefunctions. 



§ Note that there are no irreducible graphs of this type where the photon and the pion are 
attached to different nucleons, because energy and momentum conservation require an interac- 
tion between nucleons (for example, by pion exchange), and somewhere along the nucleon lines 
there is an energy flow of the order of the pion mass. The irreducible part of the diagram then 
contains the nucleon interaction, which decreases the number of separately connected pieces and 
renders the diagram higher order. At higher energies, where higher order contributions might 



become more important, it may well be necessary to account for such effects [15 
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First, the time-ordered graphs of type 2g and 2h get cancelled against recoil in the 
one-pion-exchange piece of the potential. In order to see this, consider the three diagrams 
of figure 3. These graphs all have the same spin-isospin structure: they differ only in the 
energy denominators. The first two graphs (figures 3a and 3b) arise when a diagram like 
2a is sandwiched between wavefunctions obtained from a potential whose long range part 
comes from pion exchange. They are proportional to 

1 

E(pi -q) + E{p 2 + q)- £(pi) - E{p 2 ) * 
1 1 
X E{p x -q)+w- E[p x ) + E{p 2 + q)+w-E{p 2 )_ 



w[E(p! -q) + E{p 2 + q)- E(p 1 ) - E{p 2 )} 



(15) 



w 2 V \w 

The first term corresponds to static one-pion-exchange in the potential; it is big, as 
antecipated, because these reducible diagrams have small nucleon energy denominators. 
The second term is smaller because of the additional small recoil numerator, while the 
dots sum higher orders in chiral perturbation theory. On the other hand, the remaining 
graph (figure 3c) is proportional to 

and exactly cancels recoil in the reducible diagrams. Similar cancellations can be found 
among the other time ordered diagrams. In other words, to this order in chiral pertur- 
bation theory we can omit the time ordered diagrams 2g and 2h by at the same time 
disregarding the energy dependence in the potential. 

Second, in the case of neutral pion photoproduction at threshold a number of the 
graphs in figure 2 will not contribute: those where the photon line is attached to a pion, 
and those that go like S ■ q, where q is the outgoing pion momentum. In particular, all 
the leading order graphs vanish (figure 2a), which immediately suggests that the cross 
section will be smaller than for charged pion production, and — of particular interest to 
us — more sensitive to two-nucleon contributions. Moreover, the 3-body time-ordered 
graphs (figure 2g) and the two body disconnected graphs (figure 2h) also vanish, so we 
can expect little influence of the energy- dependent part of the potential. As noted above, 
the loop graphs are finite for a neutral pion, and so there are no undetermined parameters 
to this order. 

Since all single-scattering contributions have been calculated to third order in small 
momenta, all that is left to calculate are the 3-body graphs of figure 4. In Coulomb 
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gauge (v ■ e — 0), only figures 4a and 4b survive. Hence, to $(g 3 ) in xPT, the threshold 
amplitude for neutral pion photoproduction on a nucleus is remarkably simple. We obtain 
the general formula 

•M* A | g=0 = M S 4 A + M%\ + M$ A , (17) 

where 

A^n = ^WEK (18) 

i 

M§\ = ^ 2 -J^^Y. < *a\(?® ■ & ~ tPtP)££\9 a > (19) 

l Z7r J Jn i<j Qij 

Mt - i«^E < ■ *» - rf rf^-^l^ > ■ (20) 

K Zir ) U i<j [{qij - k) 1 + M$\qij 

Sm A {k) is a generic overlap function. (See next section and appendices for details.) 

The single-scattering electric dipole amplitudes have been calculated to ?9(g 3 ) without 
an explicit isobar field, and are given by |16| 



o p _ eg A M w Ml M 2 
ho+ ~ Sti f n { m 2m* [6 + Kp) 16/?' 



E « +n ~ 8nf n [ 2m^ n 16/2 } - (21) 

Unfortunately, the single nucleon sector is not well understood. On one hand, the neutron 
amplitude has not been measured. On the other hand, the electric dipole amplitude, Eq+, 
has an interesting — and rather complicated — history, which we will briefly discuss here 
An expansion of the amplitude in powers of the pion energy was first used to derive 
a tree level "LET" |B| yielding a value E^f = -2.23 • 1(T 3 /M W +; this LET was later 



eg A r Ml Ml 



rederived |0| under an extended PCAC hypothesis, in a way that made explicit that it 
also included loop corrections, like rescattering. Experiments at Mainz [E0| and Saclay 



[ 2"T| suggested a violation of this LET. Subsequently, the data were reexamined, leading 



to the revised value Elf = (-2.0 ± 0.2) ■ 10" 3 /M W+ , a result in agreement with the 



"classical" LET |p2| . The source of the discrepancy is isospin violation; the difference of 
6.8 MeV between the pn and mr + threshold leads to a rapid variation of the amplitude 
in this region. The correct interpretation of the data depends critically on the details of 
this variation. The situation as of 1991 is reviewed in [p2[ . To further complicate the 
matter, it was then found that there are additional large finite loop contributions to E$f 



"For a review, see jL7]] . 
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at $(q 3 )] the so-called triangle graphs (see figure 2c) [[LlJ [^3|] [p^j . They are of higher 
order in the expansion in pion energy, but are proportional to M^ 1 [|T^] P5| . With these 
contributions — which clearly must be included — the i)(q 3 ) LET no longer seems to agree 
with the data. Moreover, it is clear that -Eq+ p is — at best — slowly converging. has 
now been calculated to ^(g 4 ), and shows no signs of converging. Evidently, it is difficult to 
escape the conclusion that the s-wave multipole is not a good testing ground of QCD. All 
is not lost however; recently novel p-wave LET's have been calculated, and found to have 
better convergence properties than the s- waves [ 26[| . Nevertheless, this failure of xPT in 
describing the s-waves without explicit isobar fields is a lesson that cannot be ignored 
here. Evidently the sensible thing to do is to make a best phenomenological estimate 
of the single-scattering contribution. The chiral prediction without the triangle graphs 
would appear to be a reasonable phenomenological estimate of Eq+, and so we assume 
the same for E^ n . 



V. NEUTRAL PION PHOTOPRODUCTION ON THE DEUTERON 

Here we will consider the deuteron. Deuteron phase shifts and properties are well 



described by the Bonn wave function [12]. We also give — in paranthesis — the chiral 
wave function || results. With the conventions defined in the appendices, the single- 
scattering contribution to the deuteron is given by 

E? = IImI/Z^ + E ^ S ^ k ^ = ~ 1M x W ~ 3 / M «+ (-1-38), (22) 
where Sd{k/2) is the deuteron form-factor: 

S d (k/2) = J d 3 p<f)* f (p) UP - \) = 0-722 (0.742), (23) 

evaluated with the Bonn and chiral wave functions, respectively, and we have used the 
phenomenological estimates: 

< P = "^f {— - |4(3 + «*)} = "2.24 x 10-/M.+ 

Kl n = -7?f{^W = 0.5 x 10- 3 /M n+ . (24) 



(The corresponding result from the "complete" ^(g 3 ) Eq.(21) will also be discussed 
shortly.) The 3-body contributions are readily obtained (see appendices for details). 
Figure 3a yields 
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eg A m n m 



roc JJ^ fcf 

/ — sin (^r) dr 

Jo r z 2 



87r(M w + m d )7r/|A; 
-2.20 x lQ- 3 /M w + (-2.18) 



(25) 



and figure 3b yields 
E d ^ = - 



egAm^m 



8vr(M^ + m d )27r/3 J J 



1 /X m 'r T T2r ls[n i( Z - \) kr \ 



dz dre- mr U z (- 



,1 .. r sin(,2 — i)fcr cos(z — ^)A;r 
+ w!) { 



[(z-\)krf [(z-\)krY 
= -0.43 x 10- 3 /M w+ (-0.39). 

The total is then given by 

E d = E d ss + E d ^ + E d ^ = -3.97 x lQ- 3 /M^ + 



r (z — |)A;r 
}) 



-3.95), 



to be compared to the experimental value |2jJ : 

E ex V = (_ 3i7 4 ± o.25) x 10" 3 /M U 



(26) 



(27) 



(2f 



Hence the simple picture provided by chiral symmetry does fairly well. In particular, 
the importance of the 3-body correction (charge exchange contribution) emerges as a 
consequence of chiral symmetry. It is gratifying that the results are not particularly 
sensitive to the details of the wave function. There are of course several serious caveats. 
Strictly speaking the i9(q 3 ) result without explicit isobars fairs badly, as is made clear in 
Table 1. This is, of course, a consequence of the theoretical failure in the single nucleon 
sector. Better experimental data is necessary in the single scattering sector. Currently 
new measurements of E^ are underway at Mainz and Saskatoon [11|. In principle, 



with an accurate measurement of the deuteron photoproduction amplitude one could 
extract the neutron electric dipole amplitude using our results. Of course, in order to 
be convinced of the soundness of this method one would have to calculate ^(g 4 ) 3-body 
effects in order to test the convergence properties of the nuclear matrix elements. Finally, 
we note that our final results are quite similar to results obtained some time ago, based 
on the photoproduction low-energy theorems in the impulse approximation, and assorted 
estimates of the three-body corrections [28] [29] | 30| . One might say that we have placed 
these successful results on a more sound theoretical footing by determining — by way of 
chiral power counting — the precise graphs that should dominate at threshold. 
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VI. CONCLUSION 



We have taken a pedagogical approach to pion photoproduction on nuclei in the frame- 
work of baryon chiral perturbation theory. The method presented allows one to make 
systematic use of chiral symmetry in a scattering process involving nuclei. In general, 
calculations are more involved than those of the single nucleon sector, since one must 
focus on the set of irreducible graphs, which requires use of time-ordered perturbation 
theory. In the special case of neutral pion photoproduction, the amplitude to i9(g 3 ) is 
simple, involving only tree level Feynman graphs evaluated in the heavy-fermion formal- 
ism. We evaluated the deuteron electric dipole amplitude using the Bonn and chiral wave 
functions, together with a phenomenological estimate for the single nucleon (impulse ap- 
proximation) contributions. The result of this calculation, like that of the 7r-deuteron 
scattering length ||, is in agreement with experiment. Hence the importance of 3-body 
contributions emerges in both cases as a consequence of chiral symmetry. The result is, 
of course, critically dependent on input from the single nucleon sector. Therefore, a more 
accurate determination of the nucleon electric dipole amplitudes is clearly required in or- 
der to make definite predictions for nuclei. There are many other processes which can be 
explored using this technology. Charged pion photoproduction is of great interest since 
in that case the single nucleon sector is well understood both theoretically and experi- 
mentally. Also, results — both scattering lengths and photoproduction amplitudes — for 
heavier nuclei like tritium or Helium would be a particularly novel way to explore the 
relevance of chiral symmetry in nuclei. 
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With field normalization convention, the differential cross section for 7<i — > rrd can be 
written as 



VII. ACKNOWLEDGEMENTS 



APPENDIX A: NORMALIZATION CONVENTIONS 



da 



(27T) 4 



(pi + k - p 2 - q) 



\M\ 2 d 3 p 2 d 3 q 



(A.l) 



2E X 2E 2 2u q 2u k (2vr) 3 (2tt) 



3' 



a 
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where px,P2 are the momenta of the initial and final deuterons, and q,k are the momenta 
of the outcoming neutral pion and photon, respectively. u a is the relative velocity of the 
incident particles, given by 

M « = fir- (A ' 2) 

In the center of mass frame one finds 



1 \M\ 2 

da = T^T" -W 1 ^ S(E 1 +uj k -E 2 - uj q ) q 2 dq dQ. (A.3) 

64-k 2 U a E X E 2 0J q 0J k 

The integration over q readily yields 



647T 2 E x + \k\ Ex E 2 u q uj k q (E 2 + uJ q ) 

dQ, (A.4) 



1 q \M\ 



64vr 2 \k\ (Ex + \k\) (E 2 + cu g ) 

where u k = k, and Ex, E 2 , and uo q are the energies of the deuterons and the pion, 
respectively. The slope of the differential cross section at threshold is defined as 



\k\ da, 1 \M\\ =Q 



1 dQ 647T 2 (^2 + M 2 + ( md + M J 

- 1 (A.5) 



64tt 2 {m d + M v ) 2 ' 

We can express the photoproduction amplitude in terms of rotationally invariant am- 
plitudes: 

— * — * 

M = iJ ■ a Mx + iJ • k q '■ a M 2 + iJ ■ qq- a M3 + iq-kx aM^, (A. 6) 

where a = e — (k ■ e)e. It is convenient to define an electric dipole amplitude, E d , such 
that, in Coulomb gauge, 

M\ q=0 = M d = Sn(m d + M n ) 2i(e • J) E d , (A.7) 
where J = \{a n + a p ) . It then follows that 

1*1 da _ 8 2 

7 ^ u " 3 d - ( } 
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APPENDIX B: THE S-MATRIX 

The S-matrix is defined as 

n I 

S fl = -i(2k)U\ Vi - V} ) (I] , = )M, (B.l) 

i=i J(2n)32Ei 

where n is the total number of external particles. For the 3-body process, ^yNN — >■ irNN, 
we find 

S fi = -i(2ir) A 5 A ( Pi - Pf ) MnN (B.2) 

(2tt) 9 ^2E 1 2E 2 2E[ 2E' 2 2E 1 2E n 

whereas the S-matrix for 7c? — > rrd is given by 

Sft = -t(27r)U A (p t -p f ) ; Md . (B.3) 

7 ' (2tt)« j2E d 2E' d 2E, 2E n 1 ' 



Therefore, 



Near threshold, 



E d ~ \Jm 2 d + P w m d , E' d ^m d , ^ = £ 2 = ^ = #2 ~ m > ( B - 5 ) 
and so we obtain 

-Md = ,„ ^ AW (B.6) 

APPENDIX C: FEYNMAN AMPLITUDES 

First we need the transition operator for ^NN — > irNN. Figure 3a yields 
«wv = 5 «i (-7— S£ ea 3c r c e M ) S Wl 2 _ 5„ 2 ( — e a3d r d ) B V2 + (1 <-> 2) 

Jit Q IV1 n ^Jn 

= ~ B vi S v! B vi tfi ^/2 + _ 9 l f2 2 Bv 2 B V2 e a3c^a3d + (1 <-> 2), (C.l) 

where (f = p — p* . We can make use of the relation 

B v S^B v ^2m(^a-v, U) (C.2) 



14 



to obtain 

iT NN = ^ m ?\ * • ^ ■ * ~ + (! ~ 2)- ( C - 3 ) 

Here we approximated go ~ ?o ~ an d chose Coulomb gauge ( eo = and e- k — 0). 
Using the relation between the transition operator of ^NN — > nNN and 7c? — > 7rc? , 
obtained in appendix b, we get 

(a) _ 2e^M 7r m J ■ e x 2 , , 
,Td " " (2-) 3 /. 3 (?-?)* ( 2 J - (C - 4) 

The contribution from Figure 3b is 

ilf N = B V1 (-f S£ q'y a ) B V1 % - e e a3c (q" + q'Y 

J TV " 7T 

V-M? ^ 2 ( I7^ + q)v V * tc3b Tfe2) Bv * + (! ^ 2 ) 
= f| (2^) 2 2A^(^i • if) { ^+ M l )r (-) r a V 6 2 + (1 - 2) 

" lA {2m? 2M«{J. * {P'P-tt) (C5) 



4/r y ^ v ^c^+^r 

It then follows that 

where 5" = p — f? — k. By sandwiching these transition operators between initial and final 
states, which include the effect of the nuclear wave functions, we obtain the total matrix 
element: 

M d =< f\T d \i > . (C.7) 

Finally, we obtain 



Mf = % 2_eg^n < ^ f2 _ > (ag) 



It is straightforward to check that < f 1 ■ r 2 — t^t 2 >= —2. 
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APPENDIX D: NUCLEAR MATRIX ELEMENTS 



In this appendix we obtain the coordinate space representation of our matrix elements 
and evaluate using the Bonn potential. Consider first the momentum dependent part of 

J ' ■■■ f ,q ,q j / _/ \ J ' '■ T / _l k x 



{p-p'f J ^ y ^ J (p-p>) 2 W 2' 



1 f d 3 rd 3 pd^ (J- e) e ~ l(P ~ P >r UP- £ 



47T 

(2tt 



,3 



4tt J v (2tt) 

■ k —■ ~ * 

x (J - el ?f / d*pe-^l>%(p- \)) 

= 2vr 2 J d 3 f 0}(f) (J • 6) ^ <Mr), (D.l) 
where is the spatial part of the deuteron wave function, and the identity, 

4 = J- d*r e —, (D.2) 
q 2 An J r 

has been used. The spatial part of the deuteron wave function is 

1 ,U(r) 1 „ ,~W(r). 
* W = 7C ( r + 7S Sl2<r) r »• < D ' 3 > 

where C/(r) and W(r) are the S-state and the D-state of the radial wave function, re- 
spectively, normalized such that j °° dr(U 2 + W 2 ) = 1. The spin operator is defined to 
be 

S 12 (r) = 3(^.r) (^-f)-^.^. (D.4) 

Since the S-state wave function has no angular dependence, the contribution from the 
S-state is easily evaluated: 



2 



7t 2 (J.e) / — dr / dxe~ l ~ x 
Jo r J-i 

'- 2 (/.,-, r^ sin A )dr (D .5) 



/c Jo r 2 2 
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Next, we consider the effect of the D-state part of the wave function, which consists of 
the cross term as well as the pure D-component. The cross term is given by 

J ■ € 7T r°° UW 

< J^Tpf >c ™ss = ^= J q dr— 

x JdSl e-^ r {J- eS 12 (f) + S 12 (r) J ■ e). (D.6) 

Here we can use {<7j, Oj] = 2Sij to show that 

J ■ eS 12 (f) + S 12 (f) J ■ e = 6 (e • f)(J • f) - 2(J • e), (D.7) 
and the cross term becomes 

— * 

J-e 7r f°° UW f ■% - -> -> 

< J^Tpf >Oro SS = T J q dr— J dn e-r {3 (e • f)(J ■ r) — (J ■ e)}. (D.8) 

The angular integration yields 

/ ■% - -> -> 47T 

dQ e~^ r (e-r)(J-r) = (J-e) (sin a - a cos a), (D.9) 

where a = fcr/2. Finally we have 

J-e ^ 8n 2 y/2 ^ /- 00 , C/W ro . Ax A;r A;r ,kr 2 . kr - 



< >Cross= y^ - ( 7 ' e1 X {3 sin y - 3 t cos y " ( y } sin y } - 

(D.10) 



The pure D-state contribution can be evaluated in similar fashion: 

< J^ZJy Estate =yjdr^- JdQ e^* S{ 2 (J- e) S 12 

= JL[ dr ^l [ d Q e -4-r 6 ( 2 (e . f)(J • e) - (J • e)) 
16 J r j 

(J-e) / dr— (8 sin— -4 cos— - sin- ). (D.ll) 



k 3 J r 4 2 2 2 

Evaluating the integrals with the Bonn potential yields 

— * 

J € 

< -z — =v» > s , -^ate = {J- e) 7.616 /m _l 

— * 

< >cw = (J • e) 0.051 /m- 1 

— * 
j — * 

< ^ 2 >D- State = (J • e) (-0.090) /m- 1 
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The contributions from the D-state wave function are clearly negligible relative to the 
S-state contributions. 

Finally, consider the momentum dependent part of J\A$ . Here we make use of the 
integral parametrization: 

1 r 1 dz 



f 

Jo 



{(q> - k) 2 + M 2 } q> 2 A) [{(q> - k) 2 + M 2 }z + (1 - z)q> 2 } 2 

r 1 dz 

'Jo (p + m ' 2 )2' 

where T = <f — zk, m' 2 = M 2 z(2 — z), and <f = p — ff . We can then write 

J-itf-k) e-q* f i ^ (J. /) (e ■ T) + (z - 1) ( J ■ g) g ■ I) ^ (D12) 



A;) 2 + M 2 } <f 2 Jo (I 2 + m' 2 ) 2 

The second term vanishes in the pure S-state. The Fourier transform of the first term is 

(J 2 + m' 2 ) 2 8tt 
which, in turn, yields the matrix element 

(J • I) (e ■ I) 1 f j3 -> js-, j3_, ^ ^-^ ^i(ff-tf- z %).r „-m'r a <^ ^ 

(T 2 + m' 2 ) 2 

(D.14) 

where f(r) = ^ — (J • r) (e • f) + m') . Considering the pure S-state effect, we 
obtain 



J d?re- lU - (J • f) (e • r) (J + m')} e~ m ' r , (D.13) 

matrix element 

2 4- m ) 2 o7T J 2 



{(g--£) 2 + m 2 }g> 2 ^ " k (I 2 + m' 2 ) 2 ^ 

/■l r p—m'r _ 

= 7T 2 / tte / rf 3 f f (f) (J • e) e* z -* )k * Mr) 

Jo J r 

-7T 2 f 1 dz f d 3 r<f) f (r) (J- r) (e • f) (- + m') e - m ' r e**'^ <j>i(r) 

J J T' 



jo J o r (z — gj&x 



, r £/ 2 sin [(z - i)fcr] 



-n 2 (J-e) dz dre mr U 2 (- + m'){ v 27 - v 1 2 / 12 }- 
Jo jo r [(z — ^j/crj' 3 l(z — ^)kr\ l 

(D.15) 

Finally we obtain 

< (FtMtF >5 " state =< / ' ?> 0,747 /ml (D ' 16) 

where we have used k t h = Q.685^ 1 (mass of 7r ). 
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TABLES 







n 






■i?(g 3 ) incomplete 


-2.24 


0.5 


-1.33 


-3.97 




0.96 


3.7 


3.6 


0.94 


experiment 


-2.0±0.2 (?) 


? 




-3.74±0.25 



TABLE I. The importance of the single-scattering contribution: The X-PT predictions at 
^(q 3 ) without the triangle graphs — serving as a phenomenological estimate — lead to agreement 
with the experimental value of p7f . 



21 



FIGURES 




FIG. 1. The anatomy of a matrix element. The ^'s correspond to the nuclear wave 
functions. The blob is the sum of all j4-nucleon irreducible graphs. 
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p i (a) H 2 ( b ) (c) 

FIG. 3. The energy dependence of the impulse approximation graphs (a) and (b) approxi- 



mately cancels the time-ordered 3-body diagram (c). 
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(b) 



(d) 





(e) (f) (g) (h) 

FIG. 4. Gauge invariant subset of 3-body interactions. In Coulomb gauge, only graphs (a) 



and (b) contribute to neutral pion photoproduction. 
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